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We study the collective excitations of Na 2 lr 03 in an itinerant electron approach. We consider a 
multi-orbital tight-binding model with the electron transfer between the Ir 5d states mediated via 
oxygen 2p states and the direct d-d transfer on a honeycomb lattice. The one-electron energy as 
well as the ground state energy are investigated within the Hartree-Fock approximation. When the 
direct d-d transfer is weak, we obtain nearly flat energy bands due to the formation of quasimolecular 
orbitals, and the ground state exhibits the zigzag spin order. The evaluation of the density-density 
correlation function within the random phase approximation shows that the collective excitations 
emerge as bound states. For an appropriate value of the direct d-d transfer, some of them are 
concentrated in the energy region a; < 50 meV (magnetic excitations) while the others lie in the 
energy region w > 350 meV (excitonic excitations). This behaviour is consistent with the resonant 
inelastic x-ray scattering spectra. We also show that the larger values of the direct d-d transfer are 
unfavourable in order to explain the observed aspects of Na 2 lr 03 such as the ordering pattern of the 
ground state and the excitation spectrum. These findings may indicate that the direct d-d transfer 
is suppressed by the structural distortions in the view of excitation spectroscopy, as having been 
pointed out in the ab initio calculation. 

PACS numbers: 71.10.Li, 78.70.Ck, 71.20.Be, 78.20.Bh 


I. INTRODUCTION 

The physics of 5d-based iridates has recently attracted 
much attention, since the competition between the large 
spin-orbit interaction (SOI) and the Coulomb interaction 
makes their physical properties quite different from those 
of the 3d transition metal compounds. Novel phases such 
as the topological insulator, the Weyl semimetal, and 
spin licmid have been explored extensively in these ma¬ 
terials [l|, [3- In particular, these research activities may 
have been accelerated by the waves of new discovery sup¬ 
plied by some representative materials. 

One such example is Sr 2 lr 04 , which shows the antifer¬ 
romagnetism with the spin-orbit coupled isospin j^g = 
1/2. The low-lying excitations have been detected by 
resonant inelastic x-ray scattering (RIXS), where the 
magnon peak exists at w < 0.2 eV and the exciton peaks 
emerge around w ~ 0.5 eV with substantial weights as a 
function of energy loss oj On the localized electron 

picture, the Heisenberg-type spin Hamiltonian has been 
derived by the strong-coupling theory The spin 

Hamiltonian seems to provide a good description for the 
observed magnetic excitations. Recently it has been pre¬ 
dicted that the magnon is split into two modes due to the 
interplay between Hund’s coupling and the SOI [IMl. 
Such band splitting has now been confirmed by the mag¬ 
netic critical scattering experiment [l^ and RIXS [l^ . 

As regards the itinerant electron picture, the band 
structure calculation has been carried out within the den¬ 
sity functional theory (DFT) [I3- It provides an in¬ 
sulating antiferromagnetic ground state. Recently the 
collective excitations have been investigated by intro¬ 
ducing a multi-orbital tight-binding model by present 
authors [T^, [H, [13 ■ The density-density correlation func¬ 


tion has been investigated within the random phase ap¬ 
proximation (RPA). Several bound states have emerged 
in the correlation function, which correspond well to the 
magnons and excitons in the RIXS experiment. Thus the 
weak-coupling theory based on the itinerant electron pic¬ 
ture could provide a good explanation of the excitation 
spectra, although there remains an issue whether the sys¬ 
tem really behaves like the Mott insulator or the band 
insulator 1^. 

Another fascinating example is Na 2 lr 03 , which we will 
study in this paper. It crystallizes in the space group 
C2 jm [22l - [^ , where Ir^+ ions constitute approximately 
honeycomb layer with a Na ion located at its centre as 
shown in figure [TJ It is an insulator with a temperature 
independent optical gap ^ 350 meV [13 . Although the 
exotic spin liquid ground state had been expected origi- 
nally, it is found the magnetic order sets in below Tn = 15 
K [13] . The type of the magnetic order is determined as a 
zigzag spin alignment shown in figure[2Kc) [13 [H,[H,[13 ■ 
The low-lying excitations have been detected by inelastic 
neutron scattering (INS) [11| and RIXS [11,[2§|, in which 
the magnetic and excitonic excitations are assigned. 

Several theoretical studies have already been car¬ 
ried out to explain those characteristics. The Kitaev- 
Heisenberg spin model [30l - [^ and a generic spin model 
[13 [13 have been derived by the strong coupling expan¬ 
sion based on the localized electron picture and the phase 
diagram has been examined in a wide range of parame¬ 
ter space. The zigzag alignment, unfortunately, seems 
to be realized only in rather extreme parameter values. 
Recently a spin model containing spin-spin anisotropic 
exchange couplings has been derived from the ah initio 
calculation, having led to the zigzag order in the ground 
state [13 ■ As for the itinerant electron picture, the band 
calculations based on the DFT have been carried out (l^ - 
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Figure 2: Possible magnetic orders in the honeycomb lattice; 
(a)the Neel order, (b)the stripy order, and (c)the zigzag order, 
respectively. 


It is found that quasimolecular orbitals are formed 
within a hexagon of six Ir ions, resulting in the insulat¬ 
ing ground state. Recently the self-interaction correction 
(SIC) has been done to the DFT [i^, and makes the 
zigzag spin order stabilize in the ground state. Accord¬ 
ingly the itinerant electron picture may provide a good 
starting point. 

When we turn our attention to the collective excita¬ 
tion, three characters of the intensity of the excitation 
spectrum are identified by the experiment [H, [2^ [2^ 
magnetic excitation for a; < 40 meV, excitonic excita¬ 
tion around ui ^ 450 meV, and contribution from the 
continuum part 600 meV < oj. Note that the last con¬ 
tribution is regarded as the transfer between jeff = f 
and ^ in the localized electron picture. Despite intensive 
theoretical effort, collective excitations of Na 2 lr 03 have 
not been studied yet on the itinerant electron picture. In 
this paper, we investigate the collective excitations with 
the weak-coupling theory based on the itinerant picture 
to address the excitation spectra. To this end, we utilize 
a rather simple tight-binding model instead of ab initio 
calculation to clarify the underlying mechanism of the 
low-energy excitations. Employing the Hartree-Fock ap¬ 
proximation (HFA) to the tight-binding model, we first 
calculate the one-electron energy as well as the ground- 
state energy. Then, we evaluate the density-density cor¬ 
relation function within the RPA with the help of the 
two-particle Green’s function in order to investigate the 


excitation spectra. 

We consider the electron transfer mediated through O 
2p orbitals and the direct d-d transfer between Ir ions. 
When the direct d-d transfer is weak, nearly flat and dou¬ 
bly degenerate bands come out due to the formation of 
quasimolecular orbitals as pointed out in previous stud¬ 
ies [38l - l40l| . Since one hole occupies in the t 2 g states per 
Ir ion and two Ir atoms are contained in the unit cell 
of the honeycomb lattice, the uppermost band is fully 
empty, thereby resulting in the energy gap between the 
occupied and unoccupied energy bands. The Coulomb in¬ 
teraction and magnetic order have minor influence on the 
formation of the energy gap, and hence the system may 
be called a band insulator. The zigzag magnetic order 
is found the most stable when the direct d-d transfer is 
weaker than a certain value. However, the energy differ¬ 
ences between the zigzag spin state and others such as the 
stripy, the Neel states are as small as several meV’s per 
Ir ion. Then, we calculate the density-density correlation 
function within the RPA on the zigzag order in ground 
state. In our calculation, the collective excitations come 
out as bound states as well as quasi-bound states with 
modifying the individual excitations of the electron-hole 
pair creation. With an appropriate value of the direct 
d-d transfer, we find several bound states with the exci¬ 
tation energy concentrated on w < 50 meV, which may 
be assigned as magnetic excitations. Other several bound 
states emerge in the energy region between uj ^ 350 meV 
and the bottom of the energy continuum, which may be 
assigned as excitonic excitations. These features of exci¬ 
tation spectra semi-quantitatively capture the character¬ 
istic of the observed results mentioned above. 

With increasing magnitude of the direct d-d transfer 
between Ir ions, the nearly flat bands get dispersive, but 
the system remains still insulating state even for the siz¬ 
able direct transfer. Within the HFA, the most stable 
magnetic order is, however, changed from the zigzag or¬ 
der to the Neel order when the direct d-d transfer exceeds 
a certain value. In such a situation when the direct d-d 
transfer is substantial, within the RPA on the Neel or¬ 
der, we have the bound states in the density-density cor¬ 
relation function, which are distributed in rather wide 
range of energy. Since this behavior is quite different 
from the experiment, the large direct d-d transfer is un¬ 
favourable in real materials. Actually a recent analysis of 
the tight-binding parameters by the ab initio calculation 
has claimed that the structural distortions of all types 
suppress the direct d-d transfer [s^. 

The present paper is organized as follows. In m we 
introduce a multi-orbital tight-binding model. In llHIl 
the electronic structure is evaluated within the HFA. In 
giVi we calculate the density-density correlation function 
within the RPA. The excitation spectra are evaluated at 
the special k spots. Section |3 is devoted to the conclud¬ 
ing remarks. 
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II. MODEL HAMILTONIAN 


Similar to the case of Sr 2 lr 04 , each Ir ion in Na 2 lr 03 
resides around the centre of oxygen octahedra. Due to 
the crystal electric field of IrOe, the energy level of the 
Cg orbitals of Ir atom is about 2-3 eV higher than that 
of the t 2 g orbitals. We take account of only t 2 g orbitals, 
and ignore the small lift of degeneracy arising from the 
distortion of IrOg octahedra in a first step. Thereby, 
as a minimal model, the Hamiltonian of a multi-orbital 
Hubbard model is defined on a honeycomb lattice, 


i? = ilso + +-ffkin, (2.1) 


with 

Hso = Cso E E (2.2) 

i nn'(7(T' 

Hi = U ^ ^ 

i,n 

+ E + iu' - J) 

i,n<.n'(7 

i,n^n' 

(2.3) 

Hkin = ^ ^ ( [Ti^i']n,n' d\^^di'n' <7 + H.C., (2-4) 

{i,i') n,n'<7 

where i?sOj Hii and i?kin are described by the annihila¬ 
tion {dina) and creation operators of 5d electron 

with orbital n (= yz, zx, xy) and spin a at the Ir site 
i. The {i,i') indicates the nearest neighbour sum, and 

'^ina = d^^^diyiu. 

The Hso describes the SOI of Sd electrons where L and 
S denote the orbital and spin angular momentum opera¬ 
tors, respectively. We use the value Cso = 0.4 eV in the 
following calculation. The Hi represents the Coulomb 
interaction between the Og electrons. Parameters satisfy 
U = U' + 2J [4l|. We use the values 17 = 1.4 eV, and 
J/U = 0.15 in the following calculation. These parame¬ 
ter values for Ir atom have been utilized also in Sr 2 lr 04 

[iiS. 

The iJkin stands for the kinetic energy described by 
the hopping matrix For simplicity, only transfers 

between the nearest neighbour Ir ions are taken into ac¬ 
count. There are three types of bond between the nearest 
neighbour Ir ions, and we call them as bond 1, 2, and 3 
as illustrated in figure |3] (a). Then, two kinds of electron 
transfer may contribute to Ti^i> between the adjacent Ir-Ir 
pair. The one is indirect transfer via oxygen 2p-orbitals 
and the other is direct transfer between the Ir 

5d orbitals The former could take place only be¬ 

tween the different 5d orbitals. It may be expressed in a 


matrix form with the bases n = yz, zx, xy in order: 

f 0 tp 0 

,\ tp 0 0 

\ 0 0 0 

(2.5) 

for belonging to bonds 1, 2, and 3, respectively. 

Here tp may be evaluated by 

“ ^pd-rr/ddpdi ( 2 - 6 ) 

where VpdTr{= Vnd.rr /V^) stands for the Slater-Koster mix¬ 
ing parameter between the O 2p and Ir 5d orbitals, 
and Epd denotes the charge-transfer energy from Ir 5d 
orbitals to O 2p orbitals. According to Harrison’s pro¬ 
cedure 0, Vpda and Epd may be estimated as —1.8 eV 
and 3 eV, respectively, and hence tp is estimated as 0.375 
eV. As regards the direct d-d transfer, the matrix 
may be expressed by means of the conventional Slater- 
Koster parameters Vddcr, Vdd-^ (= -‘2,Vdda/^), and VddS 
(= Vddal^^)- Notice that it has non-zero matrix elements 
between the same 5d orbitals as well as between the dif¬ 
ferent 5d orbitals. With the help of Harrison’s procedure, 
Vddrr is estimated as —0.51 eV. Since there exist several 
complications in real crystals, these values of transfer are 
regarded as providing only order of magnitude. 



rp(pd) _ 


0 0 -t 

0 0 
-tp 0 0 



III. ELECTRONIC STRUCTURE WITHIN THE 
HARTREE-FOCK APPROXIMATION 

Before going to the HFA, it is instructive to investi¬ 
gate the situation with no Coulomb interaction working. 
Without magnetic orders, we can define a unit cell con¬ 
taining Ir ions A and B shown in figure [3Ka). The cor¬ 
responding first Brillouin zone (BZ) is illustrated as a 
smaller hexagon in figure |3Kb). The one-electron energy 
is evaluated by diagonalizing 77kin + Hso for momenta in 
the BZ. 

Figure m shows the one-electron energy with momenta 
along symmetry lines. The origin of the energy is set 
at the top of the valence band. Each line is doubly de¬ 
generated. In the absence of the direct d-d transfer, the 
dispersion is flat or nearly flat for (^so = 0 or (so = 0.4 
eV, respectively (left top or bottom panel of figure H]) . 
It indicates the formation of molecular orbitals within 
the Ir hexagon. Since one hole exists per Ir site and two 
Ir ions are contained in a unit cell, the uppermost band 
should be empty, indicating a non-magnetic insulating 
ground state. Then, introduction of the direct d-d trans¬ 
fer changes the situation. As shown in the right panels 
of hgure [H it makes the bands dispersive, but the sys¬ 
tem is narrowly insulating. This situation contrasts with 
that of Sr 2 lr 04 , where the so-called “jes = 1/2” band, 
which is doubly degenerated when ignoring the Coulomb 
interaction, is half occupied, indicating that the system 
is a metal in the absence of the Coulomb interaction. 

Now we consider the situation that the Coulomb inter¬ 
action turns on, and that a certain magnetic order such 
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Figure 3: (a) Honeycomb lattice. Three types of bond repre¬ 
sented by blue, magenta, and green lines are labeled as 1, 2, 
and 3, respectively. The unit cell consists of ions A and B. In 
the magnetic ordering phases, Ir atoms A, B, C, and D form a 
unit cell, which is enclosed by a broken line, (b) Reciprocal lat¬ 
tice and corresponding Brillouin zones. The hexagon attached 
by K and M forms the first Brillouin zone in the absence of 
magnetic orders, while the rectangle enclosed by a broken 
line forms the magnetic Brillouin zone. Centre and corners of 
larger hexagon represent the reciprocal lattice points. 
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Figure 4: One-electron energy as a function of momenta along 
symmetry lines, when the Coulomb interaction is disregarded. 
The origin of energy is set at the top of the valence band. 
The conduction band is doubly degenerated. The left (right) 
panels show the results in the absence (presence) of direct d-d 
transfer. The top (bottom) panels illustrate the results in the 
absence (presence) of the SOI. 


as the Neel, the stripy, or the zigzag orders, material¬ 
izes in the ground state (see figure [2]) To describe 
these orders, four sublattices A, B, C, and D, as shown 
in figure EJa) , are introduced. The corresponding first 
magnetic Brillouin zone (MBZ) is denoted as the rect¬ 
angle enclosed by the broken line in figure EKb). With 
wave vector k in the first MBZ, we define the Fourier 


transform of annihilation operator as 

dxnaiy) = {i/N)i , (3.1) 

3 


where j runs over the lattice sites belonging to one of the 
four sublattices A, B, C, and D specified by A, and N is 
the number of Ir ions. 

With this notation, the one-electron energy Hq = 
ffkin + Hso is rewritten as 


H. = 


^4(k) i7o(k) 

k?e 




de(k), 


(3.2) 


where use has been made of abbreviations ^ = (A, n, a) 
and C — (A',n',tT'). The [iJo(k)]^,^' denotes the matrix 
element of the Fourier transform of Hq expressed in the 
^ basis. 

We carry out the HFA by following the con- 

First, 
ffi = 

1 i^3i^4)dl^dEd,^^di^, where 


ventional procedure as explained in 
we rewrite the Coulomb interaction as 


42|. 




= (n-mjCrm) with m = 1, 2, 3, and 4. By compar¬ 
ing this with (lOI) . we can determine the content of 
9{viV2'tV3Va)- Then, we replace Hi by 


3 ?iA2,?3A4 

(3.3) 

with 


= 5(?i6;C3?4) - 5(Ci6;?46); (3-4) 

where ^ = (A, i/). The {X) denotes the ground-state 
average of operator X. 

The expectation values of the electron density con¬ 
tained in has to be self-consistently determined. For 
this purpose and evaluating the single-particle energy, 
it is convenient to introduce the single-particle Green’s 
function in a matrix form with 24 x 24 dimensions, 


G(k,w)J^^^ =-\ j (T[dj(k,t)4(k,0)])e'‘"Mt, (3.5) 


where T is the time ordering operator, and X{t) = 
e'^ ‘Ae“*^ * with H' = . The Green’s function 

can be solved by diagonalizing the Hamiltonian matrix 
with 24 X 24 dimensions. Let the .^-th energy eigenvalue 
for k be i?^(k), and the corresponding wave function be 
[i7(k)]j_^. Then the Green’s function may be expressed 
as 


[G(k,aj)]^.^,=^ 

i 


[G(k)]^,4G(k)~i],.e 

w — i?^(k) ± id 


(3.6) 


where d denotes a positive convergent factor, and -|-(-) 
is taken when the energy level with if^(k) is unoccupied 
(occupied). Note that 

= 4E /[-iG(k,aj)]^.^,e‘-°"^. (3.7) 

k 
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Once we obtain the stable self-consistent solution, we 
could calculate the ground-state energy. Noting that 

^i?,(k) = (iJo)+2(i?n, (3.8) 

where the sum over (k^) is restricted within the occupied 
levels, we express the ground-state energy as 

{fl') = + (3.9) 

Here (Hq) is evaluated by using (|3.2jl and (|3.6jl : 

(^o) = ^5][17(k)-i],,e[i?o(k)]e.«[17(k)]c,^ (3.10) 

where the sum over (k£) is again restricted within the 
occupied levels. 


(a)V,,„=-0.10eV 


(b) V,,„=-0.51 eV 
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Figure 5: One-electron energy Ee(k) as a function of k along 
symmetry lines with (a) Vdda = —0.1 eV in the zigzag phase 
and (b) Vdda = —0.51 eV in the Nefl phase. U = 1.4 eV 
and U'jU = 0.7 with 2J + U' = U. Each band is doubly 
degenerated. The origin of energy is set at the top of the 
valence band. 


A. Numerical calculation 

Since the staggered mag netic moment is directing 
along the crystal a axis [^, , we assume this is the 

direction of the staggered moment for the Neel, stripy, 
and zigzag orders in the self-consistent procedure. As al¬ 
ready mentioned, the parameter values are set (^so = 0.4 
eV, 17 = 1.4 eV, and J/U = 0.15 in the following. As 
regards the transfer, we fix the strength of the indirect 
transfer by setting Vpda = —1.84 eV and Epd = 3.01 eV, 
or equivalently, tp = 0.375 eV. With evaluating (I3.10L 
we carry out the sum over k by dividing the MBZ into 
120 X 60 meshes. To achieve the convergence, the itera¬ 
tion of 100 ^ 500 times is necessary for some cases. 

For \Vdda\ < 0.23 eV, we find the zigzag order is the 
most stable one among the zigzag, Neel, and stripy or¬ 
ders, but the energy differences among those states are 
rather small. For Vdda = —0.1 eV, for instance, the en¬ 
ergy of the zigzag order is 0.007 eV per Ir ion lower than 
that of the Neel order and 0.011 eV lower than the en¬ 
ergy of stripy order. The the orbital and spin moments 
are found to be parallel to each other with (La) = ±0.265 
and (Sa) = ±0.085, hence the magnetic moment 0.435/rB- 
The latter value should be compared with the experi¬ 
mental value 0.22/rB- Figure [SKa) shows Ee(k) as a func¬ 
tion of k along symmetry lines in the zigzag order for 
Vddcj = — 0.1 eV. The dispersion curves experience weak 
dispersion, particularly for the uppermost valence band 
and the conduction band. The band gap opens with the 
size of ~ 0.5 — 0.8 eV. Note that it is similar to the one- 
electron energy shown in the left panels of figure H) The 
latter comes out from the non-magnetic state with the 
Coulomb interaction disregarded, implying that the en¬ 
ergy gap between the occupied and unoccupied states is 
originated not from the Coulomb interaction but from 
the formation of quasimolecular orbitals on Ir hexagons 
( 3 M [ 3 ^. This interpretation is different from the one 
based on the simple localized electron picture, which 


leads to the Kitaev-Heisenberg spin model in the strong 
coupling theory. Note that the small dispersions of the 
conduction band and the uppermost valence band and 
the sizable band gap are comparable with those of the 
ab-initio calculation (Figure 4 (b) in for example). 

For \Vdda\ > 0.23 eV, we find the Neel state is the 
most stable one. For, Vdda = —0.51 eV, which value may 
be estimated by the Harrison’s procedure, the orbital 
and spin moments are found parallel to each other 
with (La) = ±0.095 and {So) = ±0.012, and hence the 
magnetic moment is 0.119^b- Figure[n](b) shows Ei{k) 
as a function of k for Vdda = —0.51 eV. The dispersion 
depends considerably on k with rather small energy gap. 
Since the zigzag order is confirmed in the real Na 2 lr 03 , 
the small value of \Vdda\ is reasonable in this respect. 


IV. ELECTRON-HOLE PAIR EXCITATIONS 
WITHIN THE RPA 


To investigate the elementary excitations, we consider 
the density-density correlation function defined by 


f+-(q,a;) 
with 




(K?i«l(0]VqCe(0))e‘‘^*dt, 

(4.1) 


= (4/fV)5^4(k + q)de(k). (4.2) 


When k -f q lies outside the first MBZ, it is reduced 
back to the first MBZ by a reciprocal vector G. Since 
(q, w) is a matrix of 576 x 576 dimensions, we con¬ 
sider a representative spectral function defined by 


/(q,cc)^^[y+-(q,cc) 




(4.3) 
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To evaluate the correlation function, we introduce the particle-hole pair within the RPA. Then, the Green’s 


time-ordered Green’s function defined by 


function is expressed as 


Y^iq) = Fiq)[f + fF{q)]-^ = +f 

- si SI iSS J 

and use the fluctuation-dissipation theorem for w > 0 where I is the unit matrix, and 

Ei, 


Y+-{q) 




= —1 


Y^{q) 




Y^{q) 


where q = (q, w). The Green’s function Y^ {q) can be 
evaluated by means of the particle-hole propagator. The 
derivation is concisely summarized below and the detail 
is given in [d^ . 

We take account of the multiple scattering between 




(4.6) 


(4.7) 


- J Here the particle-hole propagator F{q) is defined as 

(4.5) 

(4-8) 


= f ^0 + w)]525JG(k, 


J 


By substituting (13.61) into the single-particle Green’s function, and carrying out the integration over ko in (14.81) , we 
have 






= + q)C/;,.(k + q)C/4,,,(k)C/4,,(k) 


[1 - ne{k + q)]n<!/(k) 


mCk -h q)[l - ni'Ck)] 


_a; — £'f(k-I-q)-I-i?^'(k)-I-i(5 w — i?£(k-|-q)-|-(k) — i(5_ 

I- 


(4.9) 


Equation dMl) contains collective modes, which ap¬ 
pear as bound states in the low-energy sector below the 
energy continuum of individual electron-hole pair exci¬ 
tations. Since F{q) is an Hermite matrix, F{q)~^ -|- T 
can be diagonalized by a unitary matrix. If an eigen¬ 
value becomes zero at w = a;B(q )5 ^^B(q) is regarded as 
the bound-state energy. Let the corresponding eigenvec¬ 
tor be Then, expanding around 

w = aJB(q), we have 


with 



[^(q)]gGi;CC' 

Lo - lob (q) + id ’ 


(4.10) 


[d^(q)]€i«;;ce 


(4.11) 






dio 




Inserting (I4.10p into the right hand side of (14.51) . we have 
the correlation function, 

y+“(g) = 27 r( 7 (q)d(a; - a;B(q)). (4.12) 


A. Numerical calculation 

We evaluate F{q) by summing over k in (14.91) with di¬ 
viding the first MBZ into 40 x 30 meshes. The bound 


states are determined by searching for lo to give zero 
eigenvalue in F{q)~^ -(-T within the accuracy of 0.001 eV. 
The corresponding intensities are evaluated from finite 
difference between uj = a;B(q) and w = a;B(q) + O.OOleV 
in (14.111) in place of the differentiation. When lo enters 
into the energy continuum of individual electron-hole pair 
excitations, we need to evaluate the imaginary part aris¬ 
ing from the denominator in (14.91) . To make a rough 
estimate, we sort each Ef(k -I- q) — Ef'(k) inside the en¬ 
ergy continuum in (14.91) into segments with the width of 
0.05 eV for 40 x 30 k-points, resulting in the histogram 
representation. Setting lo at the centre of each segment, 
we evaluate (021) and thereby dMl)- 

For \Vddcr\ < 0.23 eV, the zigzag order becomes the 
ground state as already mentioned[43. As a typical ex¬ 
ample in this region, we calculate the spectral function 
for Vddcr = —0.1 eV. Top panel in figure H] shows /(q, w) 
as a function of oo for q at T, M, and K' points. We have 
four bound states clustered in the region where lo is less 
than 50 meV for all the q-points. They may be called 
as magnetic excitations ,which correspond to the hump 
Aq in the RIXS spectra. Experimentally, INS detected a 
magnon mode below 6 meV (22| while RIXS identified an¬ 
other around 35 meV [^ . which correspond well to our 
calculated results. The lowest excitation energy at the T 
point is found around 9 meV. These magnetic excitations 
have quite a different origin from those of the Kitaev- 
Heisenberg spin model, since the former arise from the 
insulator based on the formation of quasimolecular or- 
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Figure 6: Spectral function /(q,a;) as a function of uj for q 
at r, M, and K' points; the top panel for Vdda = —0.1 eV 
and the bottom panel for Vdda = —0.51 eV. The middle panel 
is the RIXS spectra taken from for Q = (0, 0, 6.7) in the 
C2/m notation. The (5-function peaks for bound states and 
the histograms for the continuum states are convoluted with 
the Lorentzian function of the half width half maximum 0.002 
eV. The arrows indicate the lower boundary of the continuum 
spectra. 


bitals on Ir hexagons, while the latter are brought about 
from the model based on the localized electron picture in 
the strong coupling theory. There are no spectral inten¬ 
sities in the region between 0.050 eV< uj < 0.340 eV. In 
addition, we have fonnd about a dozen of bound states 
in the region between oj =0.340 eV and at the bottom 
of the energy continuum of individual electron-pair exci¬ 
tations, which may be called as the excitonic excitations 
and correspond to the RIXS spectra around the hump 
A. The continuous spectra start around w = 0.55 — 0.65 
eV, which may correspond to the RIXS spectra around 
the humps B and C. Accordingly, these capture qualita¬ 
tively the characteristic of the RIXS spectra shown in the 
middle panel of figure [Bl 

For \Vddiy\ > 0.23 eV, the Neel state becomes the 
ground state. As a typical example in this region, we cal¬ 
culate the spectra for Vdd<j = —0.51 eV, which is shown 
in the bottom panel of figure [B] The bound states are 
distributed over a wide range of energy region below the 
continuum of electron-hole pair excitations. The number 
of the bound states are 9, 4, and 6 for q at T, M, and 
K' points, respectively. The magnetic and the excitonic 
excitations are not sharply separated. Such spectral dis¬ 


tribution is quite different from the RIXS spectra shown 
on the middle panel in figure [B] 

V. CONCLUDING REMARKS 


We have studied excitation spectra in Na 2 lr 03 on the 
basis of the itinerant electron picture. We have em¬ 
ployed a multi-orbital tight-binding model with the elec¬ 
tron transfer mediated via the O 2p orbitals as well as 
the direct d-d transfer, considering only the t 2 g orbitals 
for Ir ions on a honeycomb lattice. We have calculated 
the one-electron energy as well as the ground-state en¬ 
ergy within the HFA, and then the density-density cor¬ 
relation function within the RPA. When the direct d-d 
transfer is weak, it is found that the zigzag order becomes 
the ground state, and that the energy bands have small 
dispersions, probably due to the formation of quasimolec- 
ular orbitals. The energy differences between the zigzag 
order and other orders have been, however, as small as 
several meV per Ir ion. We have obtained the collective 
excitations as bound states in the density-density cor¬ 
relation function. Magnetic excitations have been con¬ 
centrated in the energy region less than 50 meV, while 
excitonic excitations exist around lo > 350 meV, in qual¬ 
itative agreement with the RIXS spectra. When the di¬ 
rect d-d transfer has exceeded a certain value, on the 
other hand, the Neel order has become the ground state 
with the energy bands rather dispersive. The collective 
excitations have distributed over a wide range of excita¬ 
tion energy, which behavior is at variance with the RIXS 
spectra. 

The d-d transfer estimated by Harrison’s procedure is 
larger than the critical value. The above findings may 
indicate that, in Na 2 lr 03 , the direct d-d transfer is sup¬ 
pressed by the structural distortions of all types. This 
interpretation has been pointed out by a detailed anal¬ 
ysis based on the ab initio band calculation (3^ . Our 
present result lends support to this understanding from 
the point of view of the excitation spectra. To address 
this issue in more quantitative way, it may be necessary 
to calculate the excitation spectra in more realistic mod¬ 
els. Finally electron correlations beyond the HFA and 
RPA may have important effects on the excitation spec¬ 
tra. Studies along these directions are left in future. 
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